Ray chaos in optical cavities based upon standard laser mirrors 
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We present a composite optical cavity made of standard laser mirrors; the cavity consists of a 
suitable combination of stable and unstable cavities. In spite of its very open nature the composite 
cavity shows ray chaos, which may be either soft or hard, depending on the cavity configuration. 
This opens a new, convenient route for experimental studies of the quantum aspects of a chaotic 
wave field. 
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The quantum mechanics, and more generally the wave 
mechanics of systems that are classically chaotic have 
drawn much interest lately; this field is loosely indicated 
as "quantum chaos" or "wave chaos" 2.7 A, Ai ^ 
Practical experimental systems that display wave chaos 
are rare; best known is the 2D microwave stadium-like 
resonator which has developed into a very useful tool to 
study issues of wave chaos [2, El ■ Our interest is in an op- 
izcaZimplementation of all this; that would allow to study 
the quantum aspects of a chaotic wave field, such as ran- 
dom lasing, excess noise, localization and entanglement 

However, the construction of a high-quality closed res- 
onator (such as a stadium) is presently impossible in the 
optical domain due to the lack of omnidirectional mirror 
coatings with R — 100% reflectivity. The best one can do 
is to use a metal coating; however this has only R < 95% 
in the visible spectrum Of course, dielectric mul- 

tilayered mirrors can reach R = 99.999% (or more) but 
these are far from being omnidirectional. 

This leads to the consideration of an open optical cav- 
ity. One approach is to use a dielectric or semiconduc- 
tor microresonator with a deformed cross section and 
profit from (non-omnidirectional) total internal reflection 
[3, llJI ; however, such a microresonator is difficult to fab- 
ricate and control. Our approach is to construct an open 
cavity based upon standard high-reflectivity laser mir- 
rors. We will show, surprisingly, that this open cavity 
allows to generate hard chaos; a closed cavity is not re- 
quired for that. We will limit ourselves to prove that our 
system is classically chaotic; by definition, this is suffi- 
cient for a system to be wave chaotic. A proper wave- 
mechanical treatment, inclu ding the calculation of the 
spectrum, will be given later [l3l| . 

The basic idea is as follows (Fig. la). An unstable opti- 
cal cavity can be built with a concave (focussing) mir- 
ror with radius of curvature R and a convex (dispersing) 
mirror with radius of curvature r at distance I < R — r. 
This unstable cavity has exponential sensitivity to initial 
conditions [l^ but does not have mixing properties be- 
cause an escaping ray never comes back; therefore chaos 
cannot occur. We overcome this difficulty in the way il- 
lustrated in Fig. lb: a second cavity, mirror-symmetric 
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FIG. 1: (a) An unstable cavity is built of a spherical concave 
mirror (with radius R) and a cylindrical convex one (with 
radius r) at distance I < R ~ r. (b) Two unstable cavities 
are coupled to form a single cavity which is globally stable for 
L < 2R. (c) Schematic diagram of the cavity model utilized 
in the numerical analysis. For all calculations in this Letter 
we have assumed a = 0.003 and b = 0.025 where all length are 
expressed in meters fl^ . On the left mirror the position y and 
the velocity v define the geometry of the chaotic scattering. 



to the first, is utilized to recollect rays leaving the first 
cavity and eventually put them back near the starting 
point. The final design of our composite cavity is de- 
picted in Fig. Ic; L is the total length of the cavity and 
satisfies the relation L < 2R in order to assure the geo- 
metrical stability of the whole system. Depending on the 
values of R, r and I each of the two sub-cavities can be 
either stable or unstable. For I < R — r they are unstable: 
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this case is the first object of our study (see below for the 
case R—r<l<R where the two sub-cavities are stable). 
A generic ray lying in the plane of Fig. Ic and undergo- 
ing specular reflections on the cavity mirrors, will never 
leave this plane. As will be shown in this Letter, the 2D 
ray dynamics in this plane can be completely chaotic; 
consequently, from now on, we restrict our attention to 
the truly 2D cavity shown in Fig. Ic. 

The study of the chaotic properties of our composite 
cavity, starts from the analogy between geometric optics 
of a l ight ray and Hamiltonian mechanics of a point par- 
ticle [l7|. In this spirit we consider the light ray as a 
unit-mass point particle that undergoes elastic collisions 
on hard walls coincident with the surfaces of the mir- 
rors. Between two consecutive collisions, the motion of 
the point particle is determined by the free Hamiltonian 
H = p^/2 whereas at a collision the position r(i) and the 
velocity v(i) (|v(t)| = 1 m/s throughout this Letter) of 
the particle satisfy the law of reflection: 

r(t+)=r(t_), v(f+) = (l-2nn)v(t_), (1) 

where i± = i ± 0+, n is the unit vector orthogonal to 
the surface of the mirror at the point of impact and 
the second rank tensor nn has Cartesian components 
[nn]ij = nijij {i,j = 1,2). The dynamics described by 
Eqs.0J preserves both the phase space volumes and the 
symplectic property 0. The losses of a cavity due to 
finite reflectivity of the mirrors can be quantified by the 
finesse of the cavity, that is the number of bounces that 
leads to energy decay. For optical cavities realized 
with commercially available mirrors values for the finesse 
of 10^ (or even larger) can be achieved ■ In our model 
we assume a mirror reflectivity equal to 100% for all mir- 
rors and consequently we restrict ourselves to losses due 
to the finite transverse dimensions of the concave mirrors 
(parameter b in Fig. Ic). //chaos occurs, almost all tra- 
jectories will escape in the end the key question is 
whether a typical trajectory will survive sufficiently long 
that chaos is still a useful concept. 

We use the Poincare surface of section (SOS) 0| as a 
tool to display the dynamical properties of our composite 
cavity. There are several possibilities for defining a SOS; 
we choose as reference surface the left mirror, plotting y 
and Vy each time the ray is reflected by that mirror. In 
Fig. 2a we show the SOS generated by a single orbit for 
a cavity configuration such that I < R — r (geometrically 
unstable sub-cavities). Apparently, hard chaos occurs; 
the unstable periodic orbit bouncing along the Z— axis 
of the cavity is represented by an hyperbolic fixed point 
on the SOS. The explicit value of the Lyapunov expo- 
nent for this orbit can be easily calculated exploiting the 
above mentioned analogy between geometric optics and 
Hamiltonian mechanics. We recall that the magnifica- 
tion M [2 of the cavity shown in Fig. la can be easily 
calculated in terms of m, the half of trace of the ABCD 
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FIG. 2: Poincare maps obtained as described in the text. In 
both cases the cavity parameters a and b are chosen as: a — 
0.003, b = 0.025. Note that Vy ~ 0.1 corresponds to an angle 
of incidence of about 6° [T3l . (a) Hard-chaos configuration: 
R = 1, r = 0.25, I = 0.04. A single chaotic orbit is shown, 
(b) Soft-chaos configuration: R — 1, r — 0.9, I = 0.45. Three 
stable orbits (concentric ellipses) and a single chaotic one are 
shown. 



matrix of the cavity: 

M = m + {m^ -ly/^, 

where 

2(1-1) (l + ^-L 



(2) 
(3) 



Since the ABCD matrix of the cavity coincides with the 
monodromy matrix for the unstable periodic orbit 
bouncing back and forth along the Z-axis, the positive 
Lyapunov exponent Ao for such orbit is given by: 
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TABLE I: Results for Lyapunov exponents Ai for different 
cavity configurations wliicli are specified by lj{ (l^) as tlie 
length of the right (left) sub-cavity and R (r) as the radius 
of curvature of the concave (convex) mirror. All lengths are 
given in m, further details are given in the text. 



Configuration 


R 


r 


II 


Ir 


Ai 


UU 


1 


.25 


.04 


.04 


0.104 ± 0.004 


US 


1 


.90 


.05 


.30 


0.0165 ±0.0002 


SS 


1 


.90 


.45 


.45 


0.0040 ± 0.0001 


MM 


1 


.80 


.20 


.20 


0.0090 ± 0.00015 


SM 


1 


.80 


.40 


.20 


0.0096 ± 0.0001 


UM 


1 


.80 


.01 


.20 


0.0191 ±0.0004 



With the numerical values utilized for Fig. 2a we obtain 
M ~ 1.94 and Ao ~ 8.27 (in units of s^^). 

We may now ask what will happen if the two sub- 
cavities of Fig. Ic are stable (i.e. R~ r < I < R) so that 
there is no magnification. Will the chaoticity (partly) 
survive, or not? As shown in Fig. 2b we find that in 
that case the periodic orbit bouncing along the Z— axis 
is represented by an elliptical fixed point on the SOS, 
surrounded by a KAM island of stability in which three 
stable trajectories are clearly visible. Surprisingly, we 
find that despite the stability of the two sub-cavities, the 
KAM island of stability is embedded in a sea of chaotic 
trajectories with a positive Lyapunov exponent. There- 
fore depending on the values of the parameters R, r and 
I our cavity can exhibit either fully chaotic behavior, 
or soft-chaotic behavior with coexistence of ordered and 
stochastic trajectories. A rigorous theory for the calcu- 
lation of average Lyapunov exponents, entropies and es- 
cape time for open systems, has been developed by Gas- 
pard and coworkers in the last decade 0. In short, in 
chaotic open systems there exists a fractal set of never- 
escaping orbits, the so-called repeller on which quan- 
tities as the average Lyapunov exponents can be evalu- 
ated. In a Hamiltonian system with two degrees of free- 
dom the Lyapunov exponents come in pairs (Ai, A ?;), 
{i = 1,2) with Ai > A2 = and + A_^ = [ig. 
We have calculated |2fl|, |2l| the values of Ai for differ- 
ent long-living trajectories belonging to the repeller for 
different cavity configurations; the "pairs rule" (sum of 
all Lyapunov exponents equal to 0) has been confirmed 
within our numerical accuracy. The results are shown 
in Table I. Each cavity configuration is labelled as LR 
where L,R — U, M, S are labels which indicate the sta- 
bility properties of the left and right sub-cavity respec- 
tively: U = Unstable (m > 1), M = Marginally stable 
{m = 1), S ~ Stable (m < 1). In all cases we find a 
positive Lyapunov exponent which confirms that chaos 
has developed; the maximum value (Ai — 0.104) occurs 
for the UU case. We stress the fact that the values for 
the parameters characterizing the different cavity config- 
urations given in Table I are experimentally realistic (i.e. 



such mirrors are commercially available). 

Note that the value Aq — 8.27 quoted above for the 
periodic orbit bouncing along the Z— axis of the un- 
stable cavity is quite different from the average value 
Xuu ~ 0.10. We argue that this is due to the fact that 
the same unstable periodic orbit may be considered as 
an orbit of the single half-cavity (very open system: no 
chaos at all) or as an orbit of the overall composite cavity. 
Consequently the sub-cavity value for Ai, even though it 
remains obviously the same, is "diluted" into the value 
for the whole cavity. 
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FIG. 3: Escape-time as a function of the impact parameter y 
for the same hard-chaos cavity configuration utilized in Fig. 
2a, shown at different relative magnifications (Ix, 15 x and 
20 x) of the horizontal axis. The finite value (« 4000 sees) for 
each singular peak corresponds to the finite value (6 x 10'') of 
the maximum number of bounces allowed in the simulations. 
In the three pictures are clearly visible the "windows of con- 
tinuity" (arrows) and the repetition of the same pattern on 
different scales (self similarity). 

It is possible to obtain a picture of the repeller set 
by plotting the escape-time function For each tra- 
jectory starting at position y ("impact parameter") and 
horizontal velocity v on the left mirror (see Fig. Ic), we 
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calculate the time at which the last bounce occurs before 
the ray leaves the cavity. Trajectories which escape in a 
finite time (almost all), give a finite value for the escape- 
time function, whereas trapped orbits are represented as 
singular points. By definition, all initial conditions lead- 
ing to a singularity of the escape-time function belong to 
the repeller itself. In Fig. 3 we show a portion of the es- 
cape time function for a hard-chaos cavity configuration 
{UU: same parameters as in Fig. 2a). In Fig. 3a we ob- 
serve clearly three "windows of continuity" 18] (arrows) 
for which the escape time function has a small value. 
Actually, Fig. 3a as a whole results from the blow-up 
of the escape time function in a region bounded by two 
other windows of continuity which are partially visible 
on the left and the right side of the figure. Consecutive 
blow-ups of Fig. 3a are shown in Fig. 3b and Fig. 3c. 
We note that going from one picture to the next both 
the density and the height of the singular peaks increase 
thus indicating that the repeller set is dense. Moreover 
the windows of continuity act as convenient markers of 
the self-similar nature of the pattern as a whole; this self- 
similarity strongly suggests that the escape time function 
is singular on a fractal set, as expected for a repeller. The 
dense occurrence of singular points is a clear signature of 
the mixing mechanism due to the confinement generated 
by the outer concave mirrors. Note that typical escape 
times are much larger than l/Xjju, thus allowing ample 
time for chaos to develop. 

Thus, we have shown that it is possible to build, with 
commercially available optical elements, a composite op- 
tical cavity which displays classical chaotic properties. 
Despite the "local" Hamiltonian structure of its phase 
space our optical cavity is, as we had anticipated, an open 
system. Opening up a closed chaotic Hamiltonian system 
may generate transiently chaotic behavior due to the es- 
cape of almost all the trajectories and our composite 
cavity promises an easy experimental realization thereof. 
Evidence for ray chaos comes both from the computation 
of Lyapunov exponents and from the plot of the escape 
time functions. The huge density of singular points in the 
escape time functions is a strong indication that the an- 
ticipated mixing works: unstable orbits which leave one 
half-cavity are recollected by the other half-cavity until 
they become again unstable and come back to the first 
cavity. Orbits for which this process is repeated forever 
generate singularities in the escape function. 

Finally, the most convenient experimental way to re- 
alize the composite cavity as a 3D system seems to be 
by using spherical elements for focussing mirrors and a 
cylindrical element for the dispersing one. A bi-convex 
cylindrical mirror (with X as the axis of the cylinder) is 
dispersing for trajectories lying in a plane orthogonal to 
the X-axis and is neutral (flat surface) for trajectories 
in a plane containing the X-axis itself. In this case Fig. 



la represents the unstable cross section of the left sub- 
cavity. The realization of such an open chaotic cavity 
in the optical domain opens a broad perspective: many 
quantum-optics experiments can now be done on a prac- 
tical chaotic system (SulOjJ. Such experiments will greatly 
benefit from the ease of manipulation and control offered 
by the macroscopic nature of our composite cavity: work 
along these lines is in progress in our group. 

This project is part of the program of FOM and is also 
supported by the EU under the IST-ATESIT contract. 
We acknowledge J. Dingjan and T. Klaassen for stimu- 
lating discussions. 
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